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Dedicated to Michael S. Yakir 

1. Introduction 

Motivation: Problems from control theory often involve a set of physical pa- 
rameters, for instance, masses, spring constants, and damping coefficients with 
mechanical systems, or resistances, capacitances, and inductances with electrical 
circuits. The structural properties of the control system may depend crucially on 
the specific choice of concrete parameter values. In many relevant examples, a sys- 
tem is generically controllable (i.e., controllable for almost all possible parameter 
values) , but becomes uncontrollable when certain relations between the parameters 
are fulfilled. 

It has been shown for many system classes of practical interest that controlla- 
bility amounts to the torsion-freeness of a module associated to the system. For 
example, if A = K[di, . . . , d n ] for a field K and J- = C°°(R n , K), the system given 
by the linear constant-coefficient partial differential equations R{d\ 1 . . . , d n )w = 0, 
where R £ A pxq and w e T q , is controllable if and only if M = A lxq /A lxp R is 
torsion-free. Then the parametric controllability problem can be formulated as 
follows: Given R £ A pxq , where A = K(pi, . . . ,pt)[di, ■ ■ ■ , d n ] for some parameters 
px, . - .,Pt, find out whether M is generically torsion- free, and moreover, determine 
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the relations among the pi that will cause torsion elements in M. In general, we 
pose this question for a left module M over a (non-commutative) algebra A. The 
antipode notion of controllability is autonomy, which happens for systems being 
torsion. For a corresponding system module M, this means that M is annihilated 
by a non-zero ideal in A. 

To a system module M over a system algebra A one associates the transposed 
module N = N(M), defined as follows. Let the left module M be presented by 
a matrix R 6 A pxq , then N(M) is a right module with the presentation matrix 
R T eA" x P. 

Then, there is an alternative description of torsion-freeness (controllability) and 
torsion (autonomy) in the language of homological algebra. Namely, 
M is torsion-free if and only if Ext\(N(M), A) = 0, and 
M is a torsion module if and only if Ext^(M, A) = Hom^(M, A) = 0. 

For a survey of the correspondence between control systems and their system 
modules, see [5] and the references therein, in particular, the works of Oberst, 
Pommaret and Quadrat. A general approach to parametric modules including the 
case when the parameters are non-constant was introduced in jl4j [15] . The authors 
showed that so-called trees of integrability conditions, depending on parameters 
of the system, determine the control-theoretic properties of the system. These 
trees result into systems of partial differential equations and nonlinear differential 
conditions. 

The situation described above was our original motivation for studying parameter- 
dependent questions of homological algebra such as the specific problem outlined 
above. It turned out that apart from the concrete application area, it is a challeng- 
ing task for computer algebra to investigate parametric modules, and in particular, 
to get a grip on the special values of parameters that cause a qualitative change of 
structural module properties. These questions reach far beyond the limited set of 
algebras that typically arise in control theory. Roughly speaking, the problem can 
be tackled from the computational point of view for virtually every algebra that 
is accessible to Grobner basis techniques. The main idea is simple but effective: 
it consists in a careful monitoring of denominators of cofactors that appear during 
Grobner basis computation. Thus in this article we continue with the investiga- 
tions, started in the articles [9| flO] . 

Outline of the paper: In this paper, we give an algorithm for answering the 
following 

Question from Control Theory. Given a linear system S, depending on a finite 
number of parameters. Determine the control-theoretic properties (such as the de- 
composition into a controllable and an autonomous part) of S for all the values of 
involved parameters. 

Since there is, for certain system classes [2j [T01 [13] , a one-to-one correspondence 
between control-theoretic properties of a system S and the homological properties 
of an associated module M(S), we can reformulate the question as follows: 
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Question for Computer Algebra. Given a finite presentation of a parametric 
module M over a (non-commutative) algebra A, determine the properties (e.g. ho- 
mological) of M for all the values of involved parameters. 

We present detailed solutions for the bipendulum equations (Example 13. 1[) , and 
for the "two pendula mounted on a cart" problem, for both negligible friction (Sec- 
tion ^. ip and essential friction (Section l4.2p . The latter problem, to the best of our 
knowledge, has not been yet solved completely in an explicit way. We also present 
and comment on several curious examples. Namely, we show the existence of a 
non-generic controllability in the generically autonomous system (Example !3-2[) , 
and present a system, where both controllability and autonomy properties appear 
only in the non-generic situation (Example 13. 3|) . In the treatment of the case 3 
of 14.21 we illustrate the ability of our method to treat nested obstructions, that is, 
investigating sub-obstructions of a given obstruction to genericity. 

Preliminaries: Algebraically speaking, a parameter is a non-zero (and thus 
invertible) element of the ground field. In this article we deal with parameters which 
mutually commute with the elements of the algebra. In other words, the action of 
operators of the algebra on parameters is just the commutative multiplication. 

In this article, we use the following definition of a property being generic. 

Definition 1.1. Let V(pi, . . . ,p n ) be a polynomial expression in pi over some 
domain T>, on which a measure /i exists. The identity V — holds generically in 
V if 7?(6, . . . , C n ) = 0, for almost all (&,..., f„) £ V n . 

In other words, V = holds in D" \ E, where E C V n and fi(E) = 0. 

For instance, if "D — C and V vanishes on the complement of a nonzero algebraic 
set E C T>" , then V = holds generically in T> n . The Krull dimension Kr. dim 
of the coordinate ring K[V] of a variety V can be used for defining a measure on 
closed subsets E C V n by assigning fi(E) = if Kr. dimKfE 1 ] < n, and fi(E) = 1 
otherwise. 

Notations: For a matrix M, M T denotes its transposed matrix. By a{E) we de- 
note a left A-submodule, generated by the finite set F. The subscript A is dropped 
when A is commutative. For an ideal 7 in a commutative ring K[xi, . . . ,x n ], we 
denote by V(I) C K™ the set of common zeros of polynomials in I. 

2. Genericity of Grobner Bases of Parametric Modules 

Since the major role in computations (of e.g. homological properties) is played 
by Grobner bases, we investigate their behaviour in the case when a ground field 
involves parameters. 

Let IK be a field. Let A be a (non-commutative) algebra over K(pi, . . . ,p t ). 
Suppose that in this algebra the notion of algorithmic left Grobner basis exists 
(e.g. A can be a ring of solvable type |6j or, more restrictively, an Ore algebra [5]). 

Let us recall the definition of a ring of solvable type. 

Definition 2.1. Let K D K be a skew field and let R' := K \x\, . . . , x n ] be a 
commutative ring over K. Suppose that -< is a fixed term ordering on R' . Let 
R be a ring generated by {x\, ...,x n } subject to the new multiplication *. If the 
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properties 1 and 2 below hold and (R, *) is an associative ring, R is called a ring 
of solvable type. 

1. VI < i < j < n, Xi * Xj — XiXj and Xj * Xi — djXiXj + Pij, where ^ Cij e K 
and pij G i?', such that lm(py) -< a^x.,-, 

2. Vl<i<n, VaGif, d*a:i = ax^ and * a = c ai axi+p a i, where ^ c ai e if 
and Pai e K. 

Good examples of rings of solvable type are the rings of (partial) differential- 
difference operators. 

The elements of a free module A m are represented as the vectors i — X^fcLi ^i e ii 
where ti G A, and is the i-th canonical basis vector. By we denote the zero 
vector (0, . . . , 0) T G A m . The set of vectors t\, . . . , ij, for instance the set of gener- 
ators of a submodule of a free module A m , will be often identified with the matrix 
T C Mat (to x Z, A). A single vector ti corresponds to the z-th column of T and vice 
versa. 



Given a monomial well- ordering ~< on A, there are several ways to extend it to a 
monomial module ordering -<m on A m , that is, an ordering consisting of two com- 
ponents (-<, -<c), where -<c is an ordering on the components e^. In the following, 
we need a so-called term over position ordering, that is, ra^ei <m m 2 ej if and 
only if mi ~< m 2 or, if mi — to 2 , then ei -<c e j f° r monomials rrii G A. 

Recall, that a left syzygy of a finite set of elements {/i, . . . ,/ m }, fi e A, is 
a tuple (&i, . . . , b m ) T e A m , such that 61/1 + • • • + b m f m = 0. The set of all left 
syzygies of a given set of to elements is a left submodule of A m . It is often denoted 
as Syz({/i,...,/ m }). 

In this article, we work with left submodulcs, left syzygies etc. It is clearly pos- 
sible to do the same also from the right. However, two sided (bimodulc) problems 
deserve, except for the commutative fairly distinct treatment. Most (if not 

all) problems, originating from applications of e.g. control theory, are formulated 
in terms of left modules. 



2.1. Lift and Leftlnverse Algorithms. 



Proposition 2.2. Suppose that a left submodule L of a free module A m is gen- 
erated by the set of column vectors F = {/i, •••,/;} C A m . Consider the set 
F := {/i + e m+ i, ...,/; + e m+ i} and assume, that the fixed ordering -< on A m , 
naturally extended to the ordering -<i on A m+l , satisfies x a e m+ i -</ x^ej, for all 
l<i<l,l<j<m and for all a, (3. Suppose that the left Grobner basis G of F is 
finite. Then we reorder the columns of G in such a way, that the elements, whose 
first to components are zero, are moved to the left. This process is schematically 
presented in the following picture: 
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Let H = {hi, . . . ,ht} be a left Grobner basis of F. Recall that we identify F 
with the matrix (fx, ■ ■ ■ , fi) € A mxl , T with an / x t matrix over A, and H with 
the matrix (hi,---, h) eA mxt . Then 

• T is a left transformation matrix between two generating sets of F, 
that is H T = T T F T holds, 

• the columns of S form a left Grobner basis of Syz({/i, ...,/;}). 

Proof. 

i i i 

Since hj = a ik f k , we have ^ a ik (f k + e m+k ) = h i +'^2 a ik e m+k . 
fe=l fe=i k =i 

Hence, the i-th column of T is (an, . . . ,au) T , and H T — (hi, . . . ,h t ) T = T T -F T . 

Let S = {si, . . . , s r } and S :— {(0, si) T , . . . , (0, s r ) T }- Since G is a left Grobner 
basis of F, for any / in A (F) n ({0} m x A 1 ) = A (F) n e^ +1 Ae fe there exists 
g G G, such that lm(g) divides lm(/). Then lm(g) € {0}" 1 x yl', hence, by the 
property of the ordering, geGfl ({0} x A 1 ) = S. Thus S is a left Grobner basis of 
F (~i ({0}™ 1 x A') and, in particular, S generates the latter. Since 

l i i 

h(fk + e m+k ) = ^ he m +k holds if and only if ^ b k f k = 0, 

k=l fc=l fe=l 

S consists of columns (bi, . . . , h) T , which are the syzygies of the set {/i, . . . , /;}. 

□ 

Remark 2.3. Clearly, for a Noetherian algebra A the algorithm terminates. 

More generally, if the left Grobner basis of F is finite, we get the transformation 
matrix in finitely many steps. Namely, in the generalized Buchberger's algorithm for 
computing left Grobner basis, we do not consider S'-poly normals between elements 
whose leading monomials include components greater than m. 

If the algebra A is commutative, the transformation matrix property translates 
into H = F ■ T. 

We call the algorithm computing the transformation matrix as above Lift(F, H). 
Note that with this algorithm we are able to trace any computation which uses 
Grobner bases. It is worth mentioning that Proposition ^. 21 shows, that with basi- 
cally one Grobner basis computation we can get three important objects, namely 
a Grobner basis of a module, a Grobner basis of the first syzygy module and a 
transformation matrix. These three applications are sometimes called Grobner 
trinity and play a fundamental role in computer algebra. 

Many problems in control theory involve parameters, which are known to be non- 
zero, or even strictly positive, for physical reasons. However, it might happen that 
the vanishing of certain algebraic expressions in the parameters has a direct impact 
on the control-theoretic properties. Very often we observe generically controllable 
parametric systems which, for some values of parameters, become uncontrollable. 

As a further application of the algorithm Lift, we compute a left inverse of 
a given polynomial matrix in the case it exists. Below, the algorithm rmLeft- 
GroebnerBasis(M) computes the monic reduced minimal left Grobner basis of a 
submodule M, which is unique for a fixed ordering ( [6l [T] ) . 
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Algorithm 2.4. LEFTlNVERSE(matrix M) 
Input: M G Mat mx „(A) 

Output: L G Mat„ xm (A), such that L ■ M — Id nxn 
or G Mat ixi (A), if no left inverse exists 

module G := rmLeftGroebnerBasis(M) 
if G 7^ Id nxn then 

report "No left inverse exists" 

return 

endif 

module N := Transpose(M) 
module K := LlFT (N, I d n xn) 
return Transpose^) 

Proof. The algorithm LeftInverse terminates as soon as Lift does. Note that 
LM = Id n xn can happen only in the case when the monic reduced minimal left 
Grobner basis of a free submodule generated by the columns {Nj} of N = M T is 
equal to Id nxn - 

In the setup of the Lift algorithm, we use H = Id n xn- Denote by K the result of 
Lift(7V, Idnxn)- Then, by the Proposition^ Id nxn = Idj ixn = K T M T = K T N. 
Hence, for L = K T we have LM = Id n y n - O 

The existence of a left inverse (or, more generally, a generalized inverse G, such 
that G ■ M ■ G = G) often gives us the information on genericity of parameters. 
Namely, one analyzes the possible vanishing of denominators of a generalized in- 
verse, as it is done in e.g. [2]. In the special case where A = K[d] is a principal 
ideal domain, consider the module M — A lxq /A lxp R. Without loss of generality, 
we can assume R has full row rank. Then M is torsion-free if and only if there 
exists a right inverse to R. 

As we have shown, computing the inverse is a special case of computing the 
transformation matrix with the algorithm Lift. In comparison with LeftInverse, 
Lift allows us to deal effectively with more general problems. 

We call the polynomials in parameters, whose vanishing implies the failure of 
generic properties, obstructions to genericity. We can compute them as de- 
scribed above using the Lift algorithm. 

There is a need for complete information on the parametric module. It consists 
of the list of properties, computed for the generic and all the non-generic cases. 
In the context of generically controllable problems, we are interested in computing 
e.g. an annihilator of a torsion submodule for the each non-generic case. Thus, we 
need to stratify the set of obstructions. 

2.2. Stratification of Obstructions to Genericity. Let K be a field of char- 
acteristic 0. Recall that a set is called locally closed, if it is a difference of two 
closed sets. A finite union of locally closed sets is called a constructible set. 

Suppose we are given a set of polynomials P = {pi, . . . ,p n } C K.[a±, . . . , a m ], 
which are irreducible over K. 

We associate to P a set C(P) := {( = (&,, . . . , Cm) G K m | Ui=iPi(0 = °}- 
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Lemma 2.5. The set C(P) is constructible. 

Proof. Let Q := {(A', A") | A' U A" = {1, . . . , n},A' n A" = 0} be the set of all 
divisions of {1, ...,n} into two disjoint complementary subsets. Let, furthermore, 
£ :=fi\(0,{l,...,n}). Then, 

C(P)= IJ {f|VjGA , p i (O = 0, Vfc e A", p fc (0 ^ 0} = 
= IJ n(fe|jeA'}))\n(K|fceA"})) = 

(j,fe)GS 

= U n ieA 'V((Pi»\n fceA "V r «p*». 

( 3 ,fc)es 

and, indeed, we see that C(P) is a disjoint union of locally closed sets. Note that 
in C(P) there is a closed subset n<V(pi); the rest of subsets are locally closed. □ 

It is convenient to represent C(P) as a binary tree, where the vertices are the 
decision points, associated to polynomials pt, and the edges represent the logical 
conditions — 0) and (p< 7^ 0), respectively. In such a way it is easy to see, 
that starting from n elements in the set P, we will have 2™ — 1 algebraic systems 
describing the locally closed components of C (P) . 

Given two ideals I, J € K[ai, . . . , a m ], an algebraic data describing a locally 
closed set V(I) \ V(J) = V(I) \ (V(I) n V(J)) can be computed with a factorizing 
Grobner basis algorithm (e.g. [S]). Such an algorithm takes I, J as input and 
returns a list of ideals, where the zero set of the intersection of the output ideals is 
contained in the V(I) and contains the complement of the V(J) in V(I). We refer 
to this algorithm as to FactGB(7, J). 

Example 2.6. Let P = {pi,P2}, then the binary tree for C(P) consists of the 
following 3 systems of equations and inequations: {pi = 0,P2 = 0}, {pi ^ 0,p2 = 0} 
and {pi = 0,p 2 ^ 0}. 

Denote V- := V{(pi)) for i = 1,2, and V 12 := V({pi,p 2 )) = V\ n V 2 - Then, the 
decomposition of C(P) can be written as V12 W (\4 \ V2) W (V2 \ Vi) = T4 2 W (Vt \ 
V12) W (V2 \ V12), where ttJ denotes the disjoint union. 

Computationally, we need to compute the Grobner basis of an ideal J12 := 
(P11P2), and two lists Li := FACTGB({pi}, P2), obtained with the factorizing 
Grobner basis algorithm, which describe V% \ V\i. 

Given a set of polynomials {fx, . . . , f s } C K[ai, . . . , a m ], we factorize them and 
form a set of pairwisc different irreducible factors P := {pi, . . . ,p n }. We sort pi 
by using a positively graded degree ordering, starting with the smaller elements. 
With such an ordering, it is easier to compute with locally closed sets. Namely, 
the bigger elements will often reduce to simpler polynomials with respect to the 
smaller elements. Thus, also the detection of empty components (that is, systems 
with no solutions) can be achieved faster. 

Lemma 12.51 is constructive indeed. Together with the presentation of locally 
closed sets using the algorithm FactGB above, we call the whole procedure Strat- 
iFYLC(list L). It takes a finite list of irreducible polynomials on the input and 
returns a list of systems of equations and inequations, corresponding to C(P). 
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2.3. The Genericity Algorithm. Let A be a K-algebra and suppose that the 
coefficients of a given system S involve parameters pi , . . . , pt ■ We interpret the 
parameters as generators of the transcendental field extension of K and we use the 
natural K(pi, . . . ,p t )-algebra structure on A. 

Algorithm 2.7. Genericity (matrix M) 

Assume, that a monomial module ordering on the algebra A is fixed. 
Input: M e Mat m xn(A) 

Output: {hi, . . . ,h s } C K[pi, . . . ,p t ], such that if a specialization of the parameters 
implies hi(pi, . . . ,pt) = 0, then a left Grobner basis of M is different from the 
generic one 



module G := rmLeftGroebnerBasis(M) 
matrix T := Lift(M, G) 
list S, H; int i,j 
for j = 1 to I 



\\G= {.9i, ...,#}€ Mat mx£ (A) 
\\ T e Mat„ x£ (A) 



i := the leading component of 
for k = 1 to n 

if (Af tfe ^ 0) then 

if (T kj ^ 0) then S := SU DENOMINATOR^) 
end if 
end if 
end for 

end for 

if (S ^ empty list) then 

H := Factorize(S') 
H := Simplify (/f) 

else H := empty list 
return H 

Proof. The algorithm FACTORiZE(list L) returns a list of monic factors of every 
polynomial of the list L. The algorithm SiMPLiFY(list L) refines a list L by removing 
doubled appearances of same elements. We may assume it also sorts L by an 
ordering, putting with the smaller elements in the beginning of the output. 

The algorithm Genericity terminates as soon as Lift terminates. Now, we 
prove the correctness. Suppose that the leading term of cjj lies in the i-th module 
component. From the property G T — T T M T it follows, that there is a presentation 
of the element Gij G A as the sum 

n n 

G ij =J2 T kjM ik = J2Tj k Ml. 

k=l k=l 

Hence, it suffices to collect only the denominators of ^ with M ik ^ 0, since 
only such elements contribute to the leading coefficient of gj . 

If some leading coefficient of the unique generic Grobner basis vanishes for some 
specialization of parameters, then the Grobner basis under such a specialization is 
different from the generic one. □ 
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Note that with the algorithm we obtain the expressions in the parameters which 
lead to non-generic Grobner bases. In order to obtain Grobner bases under spe- 
cialization, provided by hi, one cannot use the generic Grobner basis. Instead, one 
has to compute the specialized Grobner basis from scratch. 

Suppose that the output of Genericity is the list of irreducible polynomials H. 
In practice, we exclude from H the polynomials, which do not satisfy the problem - 
specific constraints for e.g. physical admissibility like non-negativity. Then, we 
apply the algorithm STRATiFYLC(iJ) and obtain a complete stratification of a 
given system with respect to its parameters. 

2.4. Comparison with Other Methods. 

2.4.1. Comprehensive Grobner bases. Comprehensive Grobner bases (see e.g. |18j ) 
were introduced by Weispfenning and generalized to rings of solvable type by Kredel 

A comprehensive Grobner basis, by definition, is a finite subset G of a para- 
metric polynomial ideal / such that <r(G) constitutes a Grobner basis of the ideal 
generated by <r(I) under all specializations a of the parameters in arbitrary fields 

(QBI). 

The construction of a comprehensive Grobner basis follows the lines of Buch- 
berger's algorithm. However, the result will be a union of trees of ideal bases 
(called Grobner systems) , where each basis is accompanied with a set of conditions 
of parameters. Being a powerful theoretical instrument, comprehensive Grobner 
bases are quite complicated to compute. To the best of our knowledge, there is no 
implementation yet, which is able to treat serious examples. 

In our approach we separate two processes, which are unified in the comprehen- 
sive Grobner basis method. Namely, we compute the tree of sets of conditions of 
parameters after the Grobner basis and transformation matrix computations. In 
such a way we avoid repeated computations in trees of ideals and sets of conditions, 
which might occur during the computation of a comprehensive Grobner basis. 

2.4.2. The Leykin-Walther Method. The method has been formalized by Leykin 
for the case of ideals [TT] and has been generalized to modules by Walther [T7] . 
The idea behind the method has been used before, however Leykin and Walther 
formulated and proved the whole framework in a complete way. In the following, 
we reformulate the Lemma 2.3 from [17] . 

Let K be a field of characteristic 0. Given a K-algebra A, we consider parameters 
as new commutative variables and perform further computations in the K-algebra 
A := _A®kK[pi, • • ■ ,Pm\- We use in A an elimination ordering -<a for the variables 
of A. Such an ordering is characterized by the property p" 1 . . .p^" -<a t, for any 
monomial t G A and any a £ N m . 

Let G = {gi, . . . ,ge} be a reduced Grobner basis for the left submodule N C A s 
with respect the position over term ordering, induced by -<a on A s . Let, moreover, 
Qn C K\p u ...,p t ] be the ideal {p G K[pi, ...,p t ] \ pA s C N}. For 9i ^ Q N A S , 
multiply all the leading coefficients with respect to K[pi, . . . ,pt] of such and de- 
note the result by h. Let a : K[pi, ■ ■ ■ ,pt] ->Kbea specialization, then if a(h) ^ 0, 
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then <j{G) = {a(gi) 7 . . . , cr(g t )} is a Grobner basis. 

This method has some drawbacks in practice. Suppose that the number of 
parameters is big and there are many obstructions, which appear in several com- 
ponents as, say, leading coefficients by a monomial 1. This situation is typical 
for generically controllable systems. Then, using the method of Leykin-Walther, 
we are forced to compute Grobner basis of a submodule of elements as described 
above, whereas a better solution would be just to collect the leading coefficients 
in parameters. Secondly, in a similar situation we get many elements in Grobner 
basis and the analysis of the impact of obstructions, e.g. the computation of the 
stratification, becomes very involved. 

On the other hand, this method allows us to handle the cases, when the parame- 
ters satisfy algebraic identities between themselves or when there are more general 
identities, involving both variables and parameters. We believe, that this method 
will be enhanced in order to overcome the described difficulties. 



3. Implementation of Algorithms 

The described method for detecting the obstructions to genericity of paramet- 
ric modules is implemented in the procedure genericity of control theory toolbox 
CONTROL. LIB Q], which is realized as a library in the computer algebra system Sin- 
gular [4]. Singular is the specialized computer algebra system for polynomial 
computations, well-known for its high performance (especially in Grobner bases- 
related computations) and rich functionality. It uses intuitive C-like programming 
language, in which the libraries are written. It is important to mention, that Sin- 
gular is distributed under GPL license, that is, it is free for academic purposes. 

The current implementation of the procedure genericity works in a little dif- 
ferent way, compared with the Algorithm 12.71 Namely, it takes as input a matrix 
T, which is assumed to be the result of the Lift algorithm. This minor modifi- 
cation allows us to compute the data, which are independent from the choice of 
a monomial module ordering. The output of the procedure genericity is a list 
of strings and thus it is ring-independent. In the first item of the list the names 
of parameters, by which we have divided in the algorithm, are collected. Every 
further item of the list contains a single non-trivial polynomial in the parameters. 

There are several algorithms in Singular, which compute (left) Grobner bases of 
modules over commutative polynomial algebras and non-commutative Gi?-algebras 
OH] . It is recommended to use the heuristic routine groebner, which often provides 
the best match for a concrete example. For more details on Singular, consult with 
the book [5] and with the website of the system [4], which contains among other 
the online documentation. The algorithm FactGB is implemented in Singular 
and is accessible via the function f acstd. 

In the library control. lib, we have implemented several functions for sup- 
porting the research in systems and control theory. Among others, there are the 
procedures Leftlnverse and LeftKernel, their counterparts Rightlnverse and 
RightKernel, as well as canonize and iostruct. 



OBSTRUCTIONS TO GENERICITY IN STUDY OF PARAMETRIC PROBLEMS IN CONTROL THEORY 



The main purpose of the library is to provide maximal relevant information 
based on a simple input. This principle led us to the development of heuristic pro- 
cedures control and autonom, which use homological computations. Respectively, 
for systems with a full row rank presentation matrix, there are dimension-guided 
procedures controlDim and autonomDim. 

Given a system algebra and a system module over it, both procedures compute 
relevant properties of a given module from the point of view of controllability (with 
the procedure control or controlDim) or autonomy analysis (with the procedure 
autonom or autonomDim). The procedure canonize takes the output of cither 
control or autonom procedure and computes reduced and tail-reduced Grobner 
bases of the objects, thus simplifying and canonizing the output. 

We illustrate the functionality of the library and the flexibility of Singular with 
the following example. 

Example 3.1. Consider a bipendulum, that is, a system, describing a bar with 
two fixed pendula of length l\ and ti respectively (e.g. [HI [3]). The system 
algebra is a commutative algebra in variable d over a field of rationals with param- 
eters g,£i,£2, that is, Q(g, l\, l<i)\d\- A system module is presented via the matrix 



LIB "control. lib"; 
option(redSB) ; option(redTail) ; 

With the LIB command we load the library. The option commands tell SINGU- 
LAR to compute reduced bases (option(redSB) ), and also reduce not only leading 
terms, but any terms in the occurring polynomials (option(redTail) ). 

It is important to mention, that any polynomial computation in Singular re- 
quires the definition of a ground ring, 
ring rl = (0 >g) ll , 12) , (d) , (c ,dp) ; 

module RR = [d~2+g/ll, 0, -g/11] , [0, d~2+g/12, -g/12] ; 

The ring we set bears the name rl, it has Q(<7, ^i, ^2) as the ground field (0 
stands for the characteristic of a field, g, 11, 12 is a list of names for param- 
eters), and the only variable d. The last comma-separated block describes the 
monomial module ordering on rl. In this case (c,dp) means the following. We 
use the descending ordering c on the module components and the degree reverse 
lexicographical ordering dp on the monomials in the same component. 

module R = transpose (RR) ; 
list L = canonize (control (R) ) ; 
L; 

We have to transpose the module R, because Singular takes the columns of 
a given matrix presentation as the generators of a module. Here is the output of 
Singular: 

[1] : 

number of first nonzero Ext: 




. We run the following code in a Singular session. 



[2] : 



-1 



[3] : 



strongly controllable (flat) , image representation: 



[4] : 



_ [1] = [(-g*12) *d~2+(-g~2) , (-g*ll) *d~2+ (-g-2) , 
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(-11*12) *d~4+(-g*ll-g*12)*d~2+(-g~2)] 

[5]: 

left inverse to image representation: 
[6] : 

_ [1 , 1] =(-11) /(g~2*ll-g~2*12) 
_ [1 ,2] =(12)/ (g~2*ll-g~2*12) 
_[1,3]=0 
[7] : 

dimension of the system: 
[8] : 
1 

[9] : 

Parameter constellations which might lead to a non-controllable system: 
[10] : 
[1]: 
g 

[2] : 

11-12 

As one can see, in the output of the procedure we provide both textual comments 
on the properties of a system and the corresponding data. The heuristics says that 
the modules Ext^, 1 (i?, rl) of a transposed module indeed vanish for i > 1 (—1 is 
returned in this situation). Hence, the system is generically controllable (the no- 
tion of strong controllability from above coincides with classical controllability for 
systems of ordinary differential equations). Moreover, the procedure computes the 
image representation, left inverse to the image representation and the dimension of 
the system. The 10-th item is the output of the procedure genericity, that is, a 
list of strings. The polynomial obstruction to genericity in this example is t\ — £2- 
The monomial obstruction g is not physically admissible. 

Let us analyze the properties of the system in the non-generic case £\ = £2 = £■ 
We do this with the help of the following code in the same Singular session: 

ring r2 = (0,g,l) , (d) , (c,dp) ; 

module RR = [d~2+g/l, 0, -g/1] , [0, d~2+g/l, -g/1] ; 
module R = transpose (RR) ; 
list L = canonize (control (R) ) ; 
L; 

We get the following output: 
[1] : 

number of first nonzero Ext: 
[2] : 
1 

[3] : 

not controllable , image representation for controllable part: 
[4]: 

_[l] = [(g),(g),(l)*d2+(g)] 
[5] : 

kernel representation for controllable part: 
[6] : 

_[1] = [0,1] 

_[2] = [1] 
[7] : 
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obstruction to controllability 
[8] : 

_[1]-[0,1] 

_[2] = [(-l)*d2+(-g)] 
[9] : 

annihilator of torsion module (of obstruction to controllability) 
[10] : 

_[l]=(-l)*d2+(-g) 
[11] : 

dimension of the system: 
[12] : 
1 

We see that the system is not controllable, since it contains a torsion submodulc 
annihilated by {Id 2 + g). However, we give both image and kernel representations 
for the controllable part of the system and describe the obstruction to controllability 
explicitly. Now, we are interested in the autonomy analysis of this non-controllable 
system, what can be achieved with the following code: 
list A = canonize (autonom(R) ) ; 
A; 

This gives us the following output: 
[1] : 

number of first nonzero Ext: 
[2] : 


[3] : 

not autonomous 
[4] : 

kernel representation for controllable part 
[5] : 

_[1] = [0,1] 

_[2] = [(-l)*d2+(-g),-l] 
_[3] = [(g)] 
[6] : 

column rank of the matrix 
[7] : 
2 

[8] : 

dimension of the system: 
[9] : 
1 

Since the 0-th Ext module of the system module RR (in other words, Hom r2 (-R-R, r2)) 
does not vanish, the system is not autonomous. In addition, we compute a kernel 
representation for the controllable part, the column rank of the presentation matrix 
and the dimension of the system. 

Parametric systems quite often are generically controllable and contain an au- 
tonomous subsystem for some special values of parameters. In the following exam- 
ple, we show that also a generically autonomous system might be controllable in a 
non-generic case. 
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Example 3.2. Let R = K(a, b) [d] be a ring. A module N = Rj (ad + b) is generi- 
cally autonomous. However, 

if a = 0, b ^ 0, then M = and thus M is autonomous, 

if a = 0, b = 0, then M is free of rank 1 and hence M is controllable. 

A general system might specialize to controllable and autonomous system in 
non-generic cases, as the next example shows. 

Example 3.3. Let R = K(a, b)[d] be a ring. Consider a module M = R 2 / 

Generically, it is neither controllable nor autonomous, the annihilator of a torsion 

submodule is (ad + b) . 

The stratification of M with respect to parameters looks as follows: 
if a 0, b = 0, a torsion submodule of M is annihilated by (d), 

if a = 0, b ^ 0, then M is free of rank 1, 

if a = 0, b = 0, then M is free of rank 2. 

Assume, that a,b G V D K. Then the space of parameters V 2 decomposes into 
a direct sum of subspaces G W E± l+J E 2 W E 3 , where G — {(a, 6) | a 7^ 0, 6 ^ 0}, 
Si = {(a, 6) I a = 0,6 ^ 0}, E 2 = {(a, 6) | a ^ 0,_6 = 0} and E 3 = {(a,_b) \ a = 
0, b = 0}. Denote by E 1 the closure of E, then dimEi = dim £3 = 1, dim £3 = in 
T> 2 . Hence, all Ei have measure and G = T> 2 has measure 1. 

Remark 3.4. There are packages like D-modules for Macaulay2, [16] . and Ore- 
Modules for Maple, [3], which have a functionality to treat some of the problems 
above. The latter package provides the possibility to reveal dangerous parametric 
denominators via the computation of generalized inverse. 

4. Example: Two Pendula, Mounted on a Cart 

Consider the Example 5.2.28 from [12] (see also the examples and solutions to 
them in [3J) describing two pendula, mounted on a cart. 

In this example, rUj is the mass and Li is the length of the i— th pendula. Re- 
spectively, ki and di are the coefficients, characterizing the friction at the joints of 
pendula. Mq denotes the mass of the cart and g is a gravitational constant. All 
these parameters can take only non-negative values. 

Let us denote Zi := ki — rriiLig for i = 1, 2. Then the presentation matrix for a 
system module is constituted by the raws of the following matrix 

m\L\d 2 ra-iLid 2 (mi + m 2 + M )d 2 — 1\ 

mi Lld 2 + did + zi mxLxd 2 

m 2 L 2 d 2 + d 2 8 + z 2 m 2 L 2 d 2 / 

We take the transposed module of the matrix. It is convenient to consider the 
columns of the matrix above as the generators of submodule of a free module. Since 
the last generator then is just (—1, 0, 0) T , we perform reduction and simplification 
of first components with respect to this generator. In such a way we obtain much 
easier presentation matrix. 



\ 
1 ad + bj 



OBSTRUCTIONS TO GENERICITY IN STUDY OF PARAMETRIC PROBLEMS IN CONTROL THEORY 



4.1. Negligible Friction. Let us assume, that the friction is negligible (that is, 
di = and ki =0). We get the simplified presentation matrix of the transposed 
module as follows: 

/ L x d 2 -g m 1 L 1 d 2 \ 
{ L 2 d 2 - g m 2 L 2 8 2 J 
The generic reduced minimal Grobner basis is the 2x2 identity matrix. With 
the Lift algorithm we obtain the transformation matrix 



L\L 2 q2 1 miLiL 2 



-d 2 - 



g 2 Li - g 2 L 2 g g 2 m 2 L 1 - g 2 m 2 L 2 

m 2 LiL 2 ^ 2 LiL 2 ^2 _ 1 

g 2 miLi - g 2 rmL 2 g 2 L 1 - g 2 L 2 g 

L\L 2 „ 2 , L\L 2 „ L 2 

-9 H — — — -= — 



\ g 2 mi(Li-L 2 ) gm\{L\ — La) g 2 m 2 (Li - L 2 ) gm 2 (L\-L 2 ) ) 

Collecting the denominators, we can see that their lorn is m\m 2 g 2 (L\ — L 2 ). 
Since mi and g are strictly positive, the only obstruction to genericity appears 
when L\ — L 2 = 0. 

Indeed, in the case L\ = L 2 = L the generic Grobner basis is ( j^q2 m 2 

hence the system is not controllable. The torsion submodule is annihilated by the 
ideal (Ld 2 ~ g), but the system is not completely autonomous. 

4.2. Essential Friction. Now, all the parameters are strictly positive. The sim- 
plified presentation matrix of the transposed module is the following 

Lid 2 + d[d + z[ TOiiid 2 

L 2 d 2 + d' 2 d + z' 2 m 2 L 2 d 2 

where z' { := — = — % q and d', := -^j- for i = 1, 2. 

The generic reduced minimal Grobner basis is the 2x2 identity matrix. The 
output of the Algorithm 12 . 71 delivers the list of three polynomials {z[, z' 2 , P}, where 

P = L\z'\ - 2L 1 L 2 z' 1 z' 2 - L x d! x d!^ + L x d'\z[ + L 2 2 z'\ + L 2 d'\z' 2 - L 2 d! x d! 2 z' x . 

z[ — means, that ki = uiiLig. This is physically admissible situation. Let us 
analyze P for the admissibility. Indeed, P is irreducible but it has a special form, 
namely 

(1) P = (L 2 z[ - L lZ ' 2 ) 2 + (L 2 d[ - L x d! 2 ) ■ (d[z 2 - d' 2 z[). 

In particular, P vanishes if both z[ and z' 2 do, so P is admissible. 
The stratification consists of 6 cases, namely 

(1) k\ = miL\g, k 2 = m 2 L 2 g, P = 

(2) ki = miL\g, k 2 ^ m 2 L 2 g, P = 

(3) ki = miL\g, k 2 ^ m 2 L 2 g, P ^ 

(4) ki ^ miL\g, k 2 = m 2 L 2 g, P = 

(5) ki ^ miL\g, k 2 = m 2 L 2 g, P ^ 
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(6) ki ^ miLig, k 2 ^ m 2 L 2 g, P = 
The setup for Singular treatment of the cases is the following: 
LIB "control. lib" ; 

ring T = (0,g) , (ml ,m2 ,L1 ,L2 ,dl ,d2 ,kl ,k2) ,dp; 

poly P = kl~2*L2"4*m2*2-2*kl*k2*Ll"2*L2"2*ml*m2-kl*dl*d2*L2*2*m2+ 

kl*d2~2*Ll"2*ml+2*kl*g*Ll"2*L2"3*ml*m2"2-2*kl*g*Ll*L2"4*ml*m2"2+ 

k2"2*Ll"4*ml~2+k2*dl~2*L2"2*m2-k2*dl*d2*Ll"2*ml-2*k2*g*Ll~4*L2*ml~2*m2+ 

2*k2*g*Ll"3*L2"2*ml"2*m2-dl~2*g*L2"3*m2"2+dl*d2*g*Ll~2*L2*ml*m2+ 

dl*d2*g*Ll*L2"2*ml*m2-d2"2*g*Ll~3*ml"2+g"2*Ll"4*L2"2*ml"2*m2"2- 

2*g~2*Ll~3*L2~3*ml~2*m2~2+g~2*Ll~2*L2~4*ml"2*m2"2; 

poly zl = kl - ml*Ll*g; 

poly z2 = k2 - m2*L2*g; 

In particular, we can see the expression for P in terms of original variables. The 
name of a ring, where the interesting parameters live as polynomials, is T . In Sin- 
gular, we can switch between different rings and also map objects. 

Case 1). k\ = m\L\g, k 2 — m 2 L 2 g, P = 0. 
Note that these three equations describe an algebraic variety, that is a closed set. 
The Grobner basis of the ideal k\—m\L\g, k 2 — m 2 L 2 g, P is k\—m\L\g, k 2 — m 2 L 2 g, 
since P vanishes, when both k± = ui\L\g and k 2 = m 2 L 2 g. Hence, it suffices to 
plug the values for ki in the corresponding system. For this, we run the following 
code: 

ring rl = (0,g,ml,m2,Ll,L2,dl,d2,kl,k2) , (d) , (c.dp) ; 
poly zl = 0; poly z2 = 0; 
module RR = 

[ml*Ll~2*d~2+dl*d+zl, 0, ml*Ll*d~2] , 

[0, m2*L2~2*d~2+d2*d+z2, m2*L2*d~2] ; 
module R = transpose (RR) ; 
list LC = canonize (control (R) ) ; 
list LA = canonize (autonom(R) ) ; 

From the output of control and autonom procedures, we conclude, that this 
system is neither controllable nor autonomous. In particular, the torsion submod- 
ule is annihilated by (d). 

Case 2). ki = miLig, k 2 ^ m%L 2 g, P = 0. 
Here we deal with the locally closed set V{{k\ — m,\L\g,P)) \ V((k 2 — m 2 L 2 g)). 
Using the following code, we get its better description. We employ a technical trick 
by modifying a ground ring in such a way, that ki have priority over the rest of 
polynomials. In such a way during the computations the relation k\ = m\L\g will 
be used as replacing k\ with m\L\g. This is achieved by using a different ordering 
like e.g. the elimination ordering (see e.g. [5]) for ki,k 2 . 

ring T2 = (0 ,g) , (kl ,k2 ,ml ,m2 ,L1 ,L2 ,dl ,d2) , (a(l , 1) ,dp) ; 

poly zl = . . . ; poly z2 = . . . ; poly P = . . . ; // we copy them from above 
ideal 12 = P,zl; 
12 = groebner(I2) ; 
f acstd(I2,z2) ; 

The output of f acstd command gives us the only component 
[1] : 

_[l]=kl+(-g)*ml*Ll 

_[2]=k2*ml"2*Ll"4+(-g)*ml"2*m2*Ll"4*L2+m2*L2"2*dl"2-ml*Ll"2*dl*d2 
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We are able to extract e.g. k 2 from the last equation explicitly: 

, m,\L\d 2 - m, 2 L\d\ 
k 2 = m 2 L 2 g H jTI "i 

Alternatively, we can express d 2 in terms of variables nii,Li,k 2 ,di. 
Computing with substitutions, we see that this system is neither controllable nor 
autonomous. The torsion submodule is annihilated by {miL\d 2 + did). 



Case 3). k\ = miLig, k 2 ^ m 2 L 2 g, 



We use the computations of the case 2 and describe a locally closed set via the 
following system of equations and inequations 

fci = m 1 L 1 g 1 k 2 - m 2 L 2 g 0, k 2 - m 2 L 2 g / jTI «i 

In order to treat both inequations involving k 2 — m 2 L 2 g 1 we introduce a new 
parameter u (thus, u is mutually non-zero in the ground field) and plug in the 
transposed system module the fake equation k 2 — m 2 L 2 g = u. 

Also this system is generically neither controllable nor autonomous. The torsion 
submodule is annihilated by (d). Compare with the annihilator for the case 2, 
which is {miL\d 2 + did). Let us investigate, for which u the properties change. 

LIB " control. lib" ; 

ring r3 = (0 ) g,ml,m2 ) Ll,L2 ) dl,d2 ) kl,k2,u) , (d) , (c,dp) ; 
poly zl = 0; poly z2 = u; 
module RR = 

[ml*Ll~2*d~2+dl*d+zl, 0, ml*Ll*d~2] , 

[0, m2*L2~2*d~2+d2*d+z2, m2*L2*d"2] ; 
module R = transpose (RR) ; 
module S = groebner(R); 
matrix T = lift(R,S); 
genericity(T) ; 

The output of genericity delivers 

[1] : 

u,m2,L2,dl 
[2] : 

ml~2*Ll~4*u-ml*Ll~2*dl*d2+m2*L2~2*dl~2 

That is, the generic annihilator of a torsion submodule of the system sub- 
ject to constraints k\ — m\L\g — 0, k 2 — m 2 L 2 g = u ^ is indeed (d). How- 
ever, if u = k 2 — m 2 L 2 g = m i L i d ^-™z L 2 d i ^ ^ nG n-generic annihilator equals 

(miL 2 d 2 + did). This illustrates the difference between two components, corre- 
sponding to cases 2 and 3. 



Case 4). k\ ^ m\L\g, k 2 = m 2 L 2 g, P = and 



Case 5). k\ ^ m\L\g, k 2 = m 2 L 2 g, P^0. 

The simplified presentation matrix for the transposed module is symmetric, that 
is, exchanging mi <-> m 2 , Li <-> L 2 , d\ <-» d 2 and k\ <-> k 2 simultaneously does not 
change the matrix. Hence, we can take the results of case 2 respectively case 3, 
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exchange the variables and get the results for case 4 respectively case 5. 
Case 6). ki ^ miLig, k 2 ^ m 2 L 2 g, P = 0. 

Recall the special structure of a polynomial P in {T]). It is easy to see, that 
if P = and one of the two summands of P is zero, so does the other. This 
observation lead us to the first conclusion: 

Lo dU 
P = 0,if/ = 4 = -f 
L\ z[ d\ 

Going back to the original variables, it translates into 

m 2 L\ _ k 2 - m 2 L 2 g _ <h 
m\L\ k\ — ra\L\g d\ 

This is especially interesting, since the values, found in [12] for showing the non- 
generic non-controllability, were m\ = m 2 = Ma = 1, L\ = L 2 = 1, d± = d 2 = 1 
and k\ — k 2 — k. As we can see, it suffices to set mi = m%, d\ = d 2l k\ — k 2 ^ 
m 2 L 2 g and L\ — L 2 for illustrating this phenomenon. 

Let us denote by a parameter t the value of the fractions in[2l Then, 

d 2 — t ■ g?i, k 2 — t ■ ki + (m 2 L 2 — t ■ miLi)g, m 2 L\ = t ■ m\L\ 

We do the substitutions for d 2 and k 2 . As a preprocessing before Grobner 
bases, we can manipulate the generators. Consider the last generator of a trans- 
posed module, that is, the last column of the transposed presentation matrix 
(rriiLid 2 ,m 2 L 2 d 2 ) T . By multiplying the column with L 2 , we can simplify it sub- 
ject to the substitution to the column (L 2 d 2 ,tLid 2 ) T . The second generator of 
the module becomes then (0,t ■ {m\L\d 2 + did + zi)) T , from which we cancel the 
parameter t out. With the following code we perform the controllability and the 
autonomy analysis for this particular case, 
ring r6 = (0 ,g, t ,ml ,L1 ,L2,dl ,kl) , (d) , (c , dp) ; 
poly zl = kl - ml*Ll*g; 
module RR = 

[ml*Ll-2*d-2+dl*d+zl, 0, L2*d~2] , 

[0, ml*Ll~2*d~2+dl*d+zl, t*Ll*d~2] ; 
module R = transpose (RR) ; 
print (R) ; 

list LC = canonize (control (R) ) ; 
list LA = canonize (autonom(R) ) ; 

We conclude that this system is neither controllable nor autonomous. In partic- 
ular, the annihilator of the torsion submodule is the ideal (miL 2 d 2 + did + k\ — 
m\L\g). Note that in view of the equation ([2]), we obtain the equivalent symmetric 
annihilator (m 2 L 2 ,d 2 + d 2 d + k 2 — m 2 L 2 g) by e.g. multiplying the previous annihi- 
lator with the constant t. 

Now let us assume, that P = but neither of its summands vanishes. The poly- 
nomial P is quadratic with respect to any of the variables mi , m 2 ,ki,k 2l d\,d 2 and is 
quartic with respect to Li and L 2 . Let us fix one of the variables mi , m 2 , ki , k 2 , d\ , d 2 
Consider the rest of variables as parameters and compute the discriminant of a cor- 
responding quadratic equation. Since the involved variables might have only posi- 
tive real values, we obtain a condition on the discriminant of a quadratic equation. 
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If we fix mi, fci or d±, we get 4m 2 L i > k 2 — m 2 L 2 g. For fixed 7x1%, k 2 or d 2 , we obtain, 

either by a direct computation or via the symmetry, that 4m ^ L 'i > ki — m\L\g. 
Notably, both inequalities cannot become equalities simultaneously. 

d 2 

Provided 4m2 2 L a > &2 — rn 2 Lig, the explicit solution with respect to, say, d\ gives 
the following expression (recall, we use the short notation Zi — ki~ rriiLig): 

d 2 (m 2 L\z x + m\L\z 2 ) ± {m 2 L\z x - rn x L\z 2 )^J d\ - \m 2 L\z 2 

di — n 72 

2m 2 L2Z 2 

Each root corresponds to a separate system. Substituting the roots into our 
system, we obtain, that as in all previous cases, it is neither controllable nor 
autonomous. The annihilators of torsion submodules are then (2m 2 L 2 d + d 2 ± 
\J d\ — Am 2 L2Z 2 ). The annihilators with respect to mi, di,L\, z\ we obtain by the 
symmetry. 

Finally, we summarize the obtained results. 

Proposition 4.1. The complete stratification of the obstructions to genericity for 
the generically controllable system with the essential friction is obtained. All the 
components of the stratification correspond to non-controllable and non-autonomous 
systems, whose torsion submodules are annihilated by one of the ideals (for i = 1, 2) 

(miLfd 2 + did), (rriiL^d 2 + did + h - rriiLig), (d), 
and (2rriiL 2 d + d t ± y 1 d\ - 4m 4 L|(fcj - rriiLig)), provided k t < 4w J L a + rriiLig. 

5. Conclusion and Future Work 

We have investigated the parameter-dependence of structural properties (such 
as torsion-freeness) of modules over certain algebras over K(pi, . . . ,p t ), where K 
is a ground field and pi are parameters. The central idea is to keep track of all 
polynomial expressions in the pi that occur as denominators during Grobner basis 
computation. These problems have practical applications in control theory as out- 
lined in the Introduction. We have shown several nontrivial phenomena that arise 
with these questions in terms of illustrative worked examples. Our goal for the 
future is to extend this approach to the study of more general parametric module 
properties, leading to the implementation of systematic procedures for such prob- 
lems. 

In particular, one is interested in working with parameters, on which the involved 
operators act nontrivially. That is, the parameters may correspond to 
(c7-)differentiable and/or (g-)shiftable functions. Then, the field K(pi, . . . ,p t ) must 
be a differential and/or a difference field. The obstructions to genericity are then 
presented as systems of differential-difference algebraic equations (DDAE) instead 
of just algebraic equations treated in this article. Though the main principles re- 
main the same, there is a strong need for specialized techniques and systematic 
computer-algebraic support for both theoretical and implementational parts of the 
further research in this area. The case of differentiable parameters was treated in 
the articles [TU [T5] , the software package OreM odules [5] seems to be able to 
provide computational support for this case. 
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Yet another important direction of investigation is the analysis of numerical phe- 
nomena, namely inexact computations with parameters defined as floating point 
numbers or as certain inequalities. The generalization of our approach to these 
domains seems to be possible with the help of e.g. cylindrical algebraic decompo- 
sition techniques. Alternatively, one may first obtain an exact symbolic solution to 
parametric problem, say, in form of the complete stratification, and postprocess it 
with numerical or symbolical-numerical tools. 
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